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, We study the effect of long-range Coulomb interaction on the phason in spin-density wave (SDW) 

within mean field theory. In the longitudinal limit and in the absence of SDW pinning the phason 
is completely absorbed by the plasmon due to the Anderson-Higgs mechanism. In the presence of 
SDW pinning or when the wave vector q is tilted from the chain direction, though the plasmon 
' still almost exhausts the optical sum rule, another optical mode appears a,t uj < 2A{T), with small 

I optical weight. This low frequency mode below the SDW gap may be accessible to electron energy 

. loss spectroscopy (EELS). 
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I. INTRODUCTION 



I It is known for some time that the long-range Coulomb 
interaction plays a rather important role in phason dis- 
persions of both charge-density wave (CDW) and 
spin-density wave (SDW) J^. As a continuation of an 
earlier analysis of CDW [07 we study in this paper the 
■ phason spectrum in SDW in the coUisionless limit. Un- 
^ ' like in CDW [|| the phason propagator I?^ in SDW is 
O not easily accessible to neutron scattering. Further, the 
^ ' large frequency limit of in SDW does not allow to de- 
^ , fine the corresponding optical weight. Therefore in the 
present paper we analyze the dielectric function e(q, w), 
which allows us to define the optical weight. Moreover, 
not only the poles of are identical to those of D^, 
but also is directly accessible to electron energy loss 
spectroscopy (EELS). 

In sharp contrast to CDW we find, that the Anderson- 
Higgs mechanism Q operates perfectly in the longitudi- 
nal limit (i.e. when q is parallel to the chain direction) 
and in the absence of SDW pinning; the phason is com- 
pletely absorbed by the plasmon and e is identical to the 
one in the normal state. In the presence of SDW pinning 
or when q deviates from the chain direction however, an- 
other optical mode appears below the quasiparticle en- 
ergy gap representing the plasma oscillations of the nor- 
mal (uncondensed) carriers in a dielectric environment. 
The optical weight of this low frequency plasmon is al- 
ways very small, but the mode itself is well separated 
from the strong high frequency plasmon. In this gen- 
eral circumstance the Anderson-Higgs mechanism is 
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weakly broken. 

It is worth noting that we do not find any acoustic 
mode in the inverse dielectric function. In the transverse 
case (i.e. when q is perpendicular to the chain direction) 
the phason completely decouples from the density fluctu- 
ations and therefore remains acoustic, but because of the 
same decoupling none of it's features shows up in e~^. 

In Sec. II the general formalism for the dielectric func- 
tion in SDW is given. In Sec. Ill we study the pole 
structure in the longitudinal limit before moving on to 
the general case of arbitrary direction of the wavenum- 
ber in Sec. IV. Sec. V concludes with a brief summary 
of our results. 



II. DIELECTRIC CONSTANT 



We consider the strongly anisotropic Hubbard model 
as introduced by Yamaji supplemented by the long- 
range Coulomb interaction 



He = 47re^ 



^^nqn_q , (1) 



where riq is the electron density operator. Within mean 
field theory for the SDW the density-density correlation 
function including the effect of long-range Coulomb 
interaction |^ is given by 

(Kn]) = (Kn])'(l + ^(Kn])')-^ , (2) 
q- 

and 

/\r, — /\n nl^- I ([",'^A])o([i5A,rt])o 

{[n,n\) - {[n,n\)Q + U i_u([sa,5A])o 

, ^ I iar 



^Voi(-77— + 777-1— tttyI • l-i) 



Here {[n,n])' is the correlation function obtained in the 
absence of long-range Coulomb interaction, iVo is the 
electron density of states per spin at the Fermi energy, 
C = C|| + V^C±sm(p, C|| = v\\q\\, C-L = v±q±, and 
v± are Fermi velocities parallel and perpendicular to the 
chain direction respectively, and (. . .) in the second line of 
Eq.(3) means the average over (f. Finally /(C, w) is the 
generalized condensate density already defined in Ref. 
[4b]. Expressions of / in certain limits needed in the 
present paper will be given later. We note that the sec- 
ond term on the right hand side of Eq.(3) corresponds to 
the contribution of the sliding SDW condensate. 
The dielectric function defined in the usual way 

e{q,LU =l + ^{[n,n]y 4 

q 

satisfies the well known relation 
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Im(-e-i) = ^Im(Kn]) . (5) 
9 
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It also obeyes the following sum rule 

/ diVLolm{-e-^) = ^ujIIcos'' + i—f sin^ , (6) 
Jo 2 

where uj^ — 47re^iVof^ = Ane'^n/m is the plasma frqucncy 
and C0S1? = The total oscillator strength depends 

only on the direction of q. 

It is worth mentioning that the phason propagator 
in SDW is given by 

=.(2A)2{((C'-^2)/) + ^iVo(C/)2 

x[l + ^A^o(9^)]"H-' , (7) 

where {[SA,SA])c is the Coulomb-corrected correla- 
tion function of order parameter fluctuations. Comparing 
Eqs. (2-3) and (7) it is readily seen that I?^ has the same 
pole structure as {[n,n\). 

III. LONGITUDINAL LIMIT 

In the longitudinal limit C = C|| • The correlation func- 
tion ([n, n\) is much simplified and written as 

([n,n])=iVoC'(^p + C'-^')~' , (8) 

where we assume that the SDW condensate is free to 
move in the chain direction. We see that Eq. (8) has the 
single pole 

co' = u^l+e (9) 

at the plasma frquency and ([n, n]) is the same as in the 
normal state because the unpinned condensate is able to 
compensate exactly for the loss of normal carriers due to 
the SDW transition. This high frequency plasmon is un- 
damped and exhausts all the available oscillator strength 
Eq. (6). In other words the phason is completely eaten 
by the plasmon (Anderson-Higgs mechanism j^). 

This situation is slightly modified in the presence of 
SDW pinning. For simplicity let us assume that the pin- 
ning is so strong that it suppresses completely the fluctu- 
ation of the order parameter ^ A induced by the density 
fluctuation riq. This amounts to neglecting the second 
term (the contribution of the condensate) in Eq. (3). 
Then we obtain 

(Kn])„d =iVoC'(l-/)[Wp(l-/)+C'-w']"' ■ 

(10) 

Due to the imaginary part of / a weak single parti- 
cle continuum appears in Im(— e^^) for < u < |C| 
and for to > \J (2A)^ -t- with square root edge at 
LD = ^ (2A)2 -I- ^2 with a;~^ decay for lo ^ oo. The 
high frequency plasmon is modified to 
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^^ = ujl + i2AfM^)-i!^]+e , (11) 

where we used the asymptotic expansion 

/(a;»2A)«-(^f[ln(^)-i|] . (12) 

This plasmon stiffens and broadens somewhat due to the 
pinning, but remains relatively sharp and still exhausts 
almost all the oscillator strength. 

There appears however another optical mode below the 
SDW gap 2A. This low frequency mode is undamped 
because it can not decay into single particle excitations, 
and it's frequency is given by 

/M = l , (13) 

where we made use of the relation 2A <C Wp. At low 

temperature (T <C Tc) this mode is well below the gap, 
and we can use the low frequency expansion 

f{^) = h + l{^)\l-9,) + ... , (14) 

where 

fd = l-2 d(/Jcosh"^(/?[l +exp(/3Acosh(yj)]-i (15) 
Jo 

is the condensate density in the dynamic limit, 

poo 

gd = 3 d</5cosh~^ (y9[l + exp(/?Acosh(/?)]~^ , (16) 
Jo 

and /3 = 1/T. The low frequency mode is then given by 



and the corresponding optical weight is evaluated as 

2"-(2)(^)(T^ • ^'^^ 

This mode is clearly the plasma oscillation due to the 
normal carriers embedded in a dielectric with a 2A gap 
produced by the pinned SDW condensate. The optical 
weight is very small (~ (2A/a;p)^) and vanishes expo- 
nentially for T — > just like the frequency of the mode 
itself. 

If the temperature approaches Tc, the frequency of the 
lower plasmon moves close to the SDW gap, and Eq. (14) 
is solved using the limiting expression 

/(.^2A) = ;J^ii^-/. , (19) 

^ ' 2 ^l-(u;/2A)2 ^ ' 

where 

=7rA2r^(w2^A2)-3/2 (20) 
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is the condensate density in the static hmit. The fre- 
quency of the mode is given by 



a;^ = (2A)^[l-(^)2] , (21) 
and the corresponding relative optical weight is 

which vanishes like (Tc—T)^ as T approaches Tg. Clearly, 
the optical weight of this midgap plasmon is maximum 
at intermediate temperature, when the frequency of the 
mode is well separated from the gap edge, but still of the 
order of the gap value. 



IV. OPTICAL MODES FOR ARBITRARY 
DIRECTION OF THE WAVENUMBER 

When q is not parallel to the chain direction, the angu- 
lar averages in Eq. (3) are not easy to perform in general. 
However, since we want to follow the evolution of the op- 
tical modes described above in the previous Section, we 
can study the q ^ limit, keeping only the information 
about the direction of q. 

If the SDW condensate is pinned, we find that the 
results of Sec. Ill for the pinned case still apply with 
only one modification. The plasma frquency uip should 
be replaced by 

Loli^) = LolicosH + i^f sinH] , (23) 

f|l 

which determines the direction dependence of the higher 
plasma frequency. Up{'&) is still large compared to 2 A 
even in the transverse limit cosi? — 0, therefore Eqs. 
(13), (17), and (21) describing the frequency of the lower 
plasmon are valid in this case as well with no direction 
dependence. It can be seen from Eqs. (18) and (22) how- 
ever, that the optical weight of the low frequency plas- 
mon will increase by a factor of (w||/w^)^ w 10^ while -d 
increases from zero to 7r/2. Therefore the lower plasmon 
is easiest to detect in the transverse limit. 
If the SDW is unpinned, Eq. (3) reduces to 

([n,n])' = -^{C^ + Ci[l-/(a;)]} (24) 

in the long wavelength limit. The corresponding dielec- 
tric function is 

2 

s(^,Lj) = 1 - 5%{cos2 1? + (^)2 sin^ mi - /(w)]} . 

(25) 

The high frequency plasmon mode is again given by Eq. 
(23) and exhausts almost all the oscillator strength. We 
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note that although the corresponding peak in Ini(— 
is somewhat broadened in general, it becomes undamped 
in the longitudinal limit as it should according to Eq. 
(8). 

The low frequency plasmon mode is obtained from Eq. 
(25) as the solution to 

/(a;) = l + (C||/a)' . (26) 

Note, that in the transverse limit (C|| = 0) Eq. (26) is 

identical to Eq. (13) describing this mode in the pinned 
case. Otherwise the mode is always higher in the un- 
pinned case. If C|| ^ C-L and T <C Tc then the mode is 
well below the gap, and Eqs. (14-16) are used to evaluate 
it's frequency as 

a;2 = ^(2A)2[l-/, + (C||/a)2](l-3,)-i . (27) 

The relative weight of this mode is 

qri—^^r[i-u+ic\\/ami-9d)-' , (28) 

which is still small, but neither the frequency of the 
mode, nor it's weight freezes out for T ^ if C|| 7^ 0. 

In the other limit, i.e. if ^|| Ci, or T ^ Tc, we use 
Eqs. (19-20) in order to solve Eq. (26) as 



The corresponding relative optical weight is now given 

by 

Y(.:;Sn.)^t-h^(i^)[i + (Cii/a)T^ 

X [i + /s + (Cii/a)']-' , (30) 

which is small and vanishes completely in the longitu- 
dinal limit where only the high frequency plasmon is 
present as in Eq. (8). 

It is to be noted that our results for the unpinned and 
pinned cases refer to two extreme situations, namely the 
pinning strength being zero and infinite respectively. In- 
corporating finite pinning strength, or in other words dis- 
cussing the situation when the fluctuations of the order 
parameter is not completely suppressed, is beyond the 
scope of this communication. Clearly, further work in 
this direction is desirable. Howcvc!r, in the experimen- 
tally interesting transverse geometry, where the predicted 
low frequency plasmon mode is most likely to be observ- 
able, the actual pinning strength does not play a crucial 
role, since the characteristics of this mode are identical 
for the pinned and unpinned cases in the transverse limit. 



V. CONCLUDING REMARKS 

We have extended an earlier analysis of the effect of 
long-range Coulomb interaction in CDW to the phason 
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in SDW. In contrast to CDW wc find that the Coulomb 
interaction has pervasive effect on the phason. Except in 
the transverse limit where the Coulomb interaction does 
not play a role, no acoustic mode is present (Anderson- 
Higgs mechanism). This is due to the fact that in SDW 
there is no effective mass enhancement and the velocity of 
the phason would be that of the Fermi velocity. The nor- 
mal carriers, having at most the same velocity, are unable 
to screen the fast phason and only optical modes survive. 
The optical sum rule for the inverse dielectric function is 
almost completely exhausted by the high frequency plas- 
mon which is present already in the normal state. Below 
Tc however, there is a small contribution of another plas- 
mon mode at a frequency uu < 2A(T) due to the uncon- 
densed carriers in a dielectric environment, which could 
in principle be accessible to EELS, most likely close to 
the transverse limit at intermediate temperatures. 
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